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Abstract The effects of the location of the crack tip in an 
annual ring and the direction of crack propagation on the 
fracture toughness of the TR crack propagation system of 
coniferous wood (T, direction normal to the notch plane; R, 
propagation direction) were analyzed by the finite element 
method in regard of the changes in elastic modulus and 
strength within an annual ring. The critical point of the 
fracture was defined as the state where the stress of a square 
element (0.125 X 0.125 mm) in contact with the crack 
tip equals the tensile strength. The distribution of specific 
gravity was measured by soft X-ray densitometry. The 
elastic moduli in the T and R directions were estimated 
by the sound velocity. The tensile strengths in the T and 
R directions were measured by the tensile test using 
small specimens of 1mm span length. Regarding the vari¬ 
ability of fracture toughness (X IC ), the experimental and 
calculated results had the same tendency. Therefore, it was 
concluded that the variability of K lc is caused by the (1) 
heterogeneity of the elastic modulus and strength within an 
annual ring; and (2) changes in the degree of stress concen¬ 
tration at the crack tips, according to the direction of crack 
propagation. 
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Introduction 

In previous reports 1,2 the variability of fracture toughness 
(K 1C ) in an annual ring was measured in the TR crack 
propagation system (when T refers to the direction normal 
to the notch plane and R to the propagation direction) of 
sugi (Cryptomeria japonica ) and spruce (Picea sitchensis). 
As a result, the critical stress intensity factors (K lc ), which 
were calculated from the load at the initiation of a crack, 
varied according to the location of the crack tip in a single 
annual ring and the direction of crack propagation 
(i.e., from pith to bark or from bark to pith). When the 
crack tip was located in earlywood, the K 1C value increased 
with the decrease in distance between the crack tip and the 
annual ring boundary in front of it; and when the crack tip 
was located in latewood, the K 1C was relatively small. 
It is believed that the variability of K lc by the crack 
tip position in an annual ring is due to: (1) the difference in 
cell shape or cell-wall thickness around the crack tip; 
and (2) the difference in stress concentration caused by 
the location of the crack tip and the direction of crack 
propagation. 

The finite element method (FEM) is useful for calculat¬ 
ing stress distributions near the crack tip. In the field of 
fracture mechanics of wood science, many finite element 
analyses have been carried out, and studies on crack initia¬ 
tion or the propagation mechanism have been reported. 3 ” 8 
Flowever, there is little research that takes into consider¬ 
ation the heterogeneity of wood (e.g., unevenness of me¬ 
chanical properties in an annual ring) except studies on the 
stress distribution in cell-wall substances 9,10 or around the 
knots. 4,7,8 

In this report, the variations of K IC according to 
the location of the crack tip in an annual ring and the 
direction of crack propagation were calculated by the FEM, 
taking into account the heterogeneity of cell properties. 
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The findings were then compared with the results of 
experiments. 


Materials and methods 

Specimens 

Spruce (Picea sitchensis Carr.), sugi (Cryptomeria japonica 
D. Don), and akamatsu (Pinus densiflora Sieb. et Zucc.) 
were used in the experiments and the FEM calculations. 
Specific gravity, moisture content, and annual ring width of 
the specimens are shown in Table 1. The single-edge 
notched-type tensile specimen of the TR crack propagation 
system was employed. Specimens were designed so all 
notch tips exist in the same annual ring. Figure 1 illustrates 
the shapes and sizes of two kinds of specimen: pith-side 
notched type and bark-side notched type. The specimen 
width was 25 mm, its thickness 10 mm, and the notch length 
11mm. Using a thin circular saw, the notch was made 
initially to 10 mm depth and was then extended 1mm more 
by a razor blade. The position of the notch tip in the annual 
ring was determined by microscopy just before the fracture 
test. 


Fracture toughness test 

The fracture toughness tests were performed at 20°C, 65% 
relative humidity using an electrohydro servo testing ma¬ 
chine (Servopulser EF1F-EA5; Shimadzu, Kyoto, Japan). 
Specimens were held by a pin-loading system. The cross¬ 
head speed was 1 mm/min. The crack opening displacement 
(COD) was measured using a clip gauge. 

The fracture toughness, (K 1C ), was calculated by the fol¬ 
lowing isotropic equation. 11 

K lc = o c -JnaF(a/W), (1) 

F(a/W) = 1.12 - 0.231 (a/W) + 10.55 (a/wf 
- 21.72 (a/W) 3 + 30.39 (a/W) 4 

where a c is the nominal critical stress, a is the notch length, 
W is the specimen width, and F(alW) is a shape factor. In 
the present work, the proportional-limit stress was o c , be¬ 
cause we recognized visually that the load of crack initiation 
agreed with the proportional-limit load in the load-COD 
diagrams. 

Distribution of specific gravity and mechanical properties 
in a single annual ring 

The distribution of specific gravity in an annual ring was 
measured with soft X-ray densitometry. Softex Type K and 
Joyce Loebl Microdensitometer 3CS were used as the soft 
X-ray irradiating device and densitometer, respectively. 
The results are shown in Fig. 2. The unit of the horizontal 
axis is defined as shown in Fig. 3. 


Table 1. Specific gravity, moisture content, and investigated annual 
ring width __ 


Species 

Specific 

gravity 

Moisture 
content (%) 

Annual ring 
width (mm) 

Sugi (Cryptomeria 
japonica D. Don) 

0.36 (0.00) 

12.4 (0.1) 

5.0 (0.3) 

Spruce (Picea 
sitchensis Carr.) 

0.44 (0.00) 

12.1 (0.1) 

4.1 (0.2) 

Akamatsu (Pinus densiflora 
Sieb. et Zucc.) 

0.55 (0.00) 

12.5 (0.0) 

4.8 (0.3) 


Numbers in parentheses are standard deviations 




Fig. 1 . Test specimens. Left Pith-side notched. Right Bark-side 
notched. Units are millimeters 


The ultrasonic pulse velocity propagating in the material 
was measured, and the distributions of elastic moduli in a 
tangential direction (£ T ) and a radial direction (E R ) within 
an annual ring were estimated by Eq. (3). 12 

E = qv 2 (3) 

where v is the sound velocity; E is the elastic modulus; and 
q is the density. The specimen sizes were 1 (T) X 1 (R) X 10 
(L) mm [for the measurement of E T (T-specimen)] and 1 
(R) X 5 (T) X 5 (L) mm (for the measurement of E R (R- 
specimen)]. Both the T and R specimens were continuously 
obtained from the same single annual ring, which was inves¬ 
tigated by the fracture toughness test. A supersonic detec¬ 
tor (Ultrasound Tester SM95; Tokyo Keiki, Tokyo, Japan) 
was used for measuring the sound velocity. The diameter of 
the transducer was 14mm. Figure 4 shows the relations 
between specific gravity and the elastic moduli (E T and E R ) 
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Relative position in an annual ring (%) 


Fig. 2. Distributions of the specific gravity in a 



single annual ring 



Relative position in an annual ring (%) 


Fig. 3. Definition of the relative 
position of the crack tip in an 
annual ring in this study. A, 
annual-ring width; B, radial dis¬ 
tance from the ring boundary of 
the previous year 



obtained by the sound velocity method. Because the elastic 
modulus of wood is generally expressed as the exponential 
power of specific gravity, 13 the distributions of E T and E R in 
an annual ring were obtained using the results of Fig. 2 and 
the regression equations shown in Fig. 4. The results are 
shown in Fig. 5. 

The distribution of elastic modulus in a longitudinal di¬ 
rection (E l ) within an annual ring was measured by the 
oscillating lead method. The specimens were 0.2 (R) X 10 
(T) X 50 (L) mm. 

The distribution of tensile strength in a tangential 
direction (c tX ) and a radial direction (a tR ) within an annual 
ring were obtained from the results of an actual tensile test. 
The specimens were the same as those used for measuring 
sound velocity. Reinforcing metal plates were bonded as 
the gripping parts of the specimens. The tensile testing 
was carried out at a crosshead speed of 0.5mm/min. Figure 
6 shows the relations between the specific gravity and 
tensile strength (cr tT and a tR ) of sugi, spruce, and akamatsu. 
The strength of wood is also expressed in the exponential 
power of the specific gravity. 14 The distributions of cr, x 
and er tR in an annual ring are shown in Fig. 7 using the 
results of Fig. 2 and the regression equations shown in 
Fig. 6. 


FEM calculation of fracture toughness 

The FEM was applied to analyze how the fracture tough¬ 
ness is affected by the difference in stress distribution that 
results from the location of the crack tip in a single annual 


ring and the direction of crack propagation. The program 
used was a universal FEM program (HITAC ISAS II), 
which is a library program of the Computer Center of the 
University of Tokyo. The analyses were carried out with a 
super computer (HITAC S-820). 

Figure 8 shows the finite element meshes used for the 
numerical calculations. The shaded portion in Fig. 8 corre¬ 
sponds to the analyzed finite element meshes. There were 
4534 elements and 4483 grids. The smallest elements (i.e., 
the elements near the crack tip) are 0.125 X 0.125 mm 
squares. The displacement of segment AB in Fig. 8 was 
fixed in the T direction, and forced displacement was pro¬ 
vided to the line DE in Fig. 8 in the T direction, realizing the 
condition of crack opening. 

The annual ring width was 5 mm, divided into 10 sections 
with different properties (elastic modulus and shearing elas¬ 
tic modulus) from each other. By shifting the property 
groups in the R direction step by step, we attempted to 
reproduce the differences in the location of the crack tips 
within an annual ring. The effect of the curvature of the 
annual ring was not considered in these FEM meshes. The 
elastic modulus (E T and E R ) and tensile strength (cr tT and 
a tR ) of each section in an annual ring were obtained from 
the above-mentioned procedures. The shearing elastic 
modulus (Ct-tr) was derived from Eq. (4). 13 

G tr = 0.003£ l (4) 

It is assumed that Poisson’s ratios (v RT , v LX , and v LR ) are 
constant, and the values noted by Fushitani et al. 13 were 
used. 
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Fig. 4. Relations between modu¬ 
lus of elasticity and specific 
gravity 




Specific gravity Specific gravity 
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Fig. 5. Distributions of the predicted modulus of elasticity in a single annual ring. Filled circles, E T ; open circles, E R 


Table 2. Mechanical properties in a single annual ring used for finite 
element method analysis of sugi wood 


Relative position 
in an annual ring 
(%) 

(X10 3 kgf/cm 2 ) 



(kgf/cm 2 ) 

£ t 

e r 

E h 


(7,t 

^tR 

0-10 

0.966 

2.20 

53.1 

0.159 

23.5 

32.2 

10-20 

0.982 

2.23 

53.6 

0.161 

23.8 

32.5 

20-30 

0.958 

2.18 

52.9 

0.159 

23.3 

32.0 

30-40 

1.13 

2.55 

57.7 

0.173 

27.1 

35.2 

40-50 

1.20 

2.69 

59.5 

0.179 

28.6 

36.4 

50-60 

1.36 

3.03 

63.6 

0.191 

32.1 

39.2 

60-70 

1.74 

3.81 

72.2 

0.217 

40.1 

45.1 

70-80 

2.78 

5.92 

92.2 

0.277 

61.5 

59.0 

80-90 

3.09 

6.54 

97.4 

0.292 

67.8 

62.7 

90-100 

6.02 

12.2 

138 

0.414 

124 

91.8 


Table 4. Mechanical properties in a single annual ring used for finite 
element method analysis of akamatsu wood 


Relative position 
in an annual ring 
(%) 

(X10 3 kgf/cm 2 ) 



(kgf/cm 2 ) 

E t 

e r 

e l 

Gtr 

0.x 

°tR 

0-10 

2.04 

3.85 

74.6 

0.224 

48.1 

56.8 

10-20 

2.23 

4.12 

78.6 

0.236 

51.1 

60.4 

20-30 

2.26 

4.16 

79.5 

0.239 

51.6 

60.9 

30-40 

2.05 

3.86 

74.8 

0.224 

48.3 

57.0 

40-50 

2.15 

4.01 

77.2 

0.232 

50.0 

59.0 

50-60 

2.30 

4.22 

80.5 

0.241 

52.3 

61.7 

60-70 

2.82 

4.95 

91.6 

0.275 

60.4 

71.2 

70-80 

7.95 

11.1 

177 

0.530 

125 

147 

80-90 

18.7 

21.7 

304 

0.911 

228 

266 

90-100 

21.0 

23.7 

327 

0.98 

248 

288 


E T , E r , E l , modulus of elasticity in tangential, radial, and longitudinal 
directions; ^TRj modulus of shearing elasticity; a tT , cr, R , tensile strength 
in tangential and radial directions; v LR , v LT , v RT , Poisson’s ratio 
v LR , 0.40; v LT , 0.60; v RT , 0.90 


E t , £ R , E l , modulus of elasticity in tangential, radial, and longitudinal 
directions; G TR , modulus of shearing elasticity; cr tT , ct, r , tensile strength 
in tangential and radial directions; v LR , v LT , v RT , Poisson’s ratio 
r'LR, 0.40; v LT , 0.60; v R -c, 0.65 


Table 3. Mechanical properties in a single annual ring used for finite 
element method analysis of spruce wood 


Relative position 
in an annual ring 
(%) 

(Xl0 3 kgf/cm 2 ) 



(kgf/cm 2 ) 


Ey 

e r 

e l 

(?TR 

OtT 

^tR 

0-10 

2.00 

2.38 

67.2 

0.202 

52.0 

59.5 

10-20 

2.00 

2.38 

67.2 

0.202 

52.0 

59.5 

20-30 

2.10 

2.46 

69.4 

0.208 

53.9 

61.2 

30-40 

2.31 

2.60 

73.3 

0.220 

57.2 

64.4 

40-50 

2.48 

2.72 

76.5 

0.229 

59.9 

67.0 

50-60 

2.61 

2.81 

78.9 

0.237 

62.0 

69.0 

60-70 

3.07 

3.10 

86.8 

0.260 

68.9 

75.4 

70-80 

3.75 

3.51 

97.8 

0.293 

78.5 

84.3 

80-90 

4.89 

4.13 

115 

0.344 

93.4 

97.6 

90-100 

7.55 

5.39 

148 

0.445 

124 

124 


Er, £ r , £ l , modulus of elasticity in tangential, radial, and longitudinal 
directions; G TR , modulus of shearing elasticity; a tT , a®, tensile strength 
in tangential and radial directions: v, R , v, T . v RT , Poisson’s ratio 
v LR , 0.37; v LT , 0.47; v RT , 0.44 


The critical point of fracture was assumed to be the load 
where the stress of one element in contact with the crack tip 
(Fig. 8B) equaled the tensile strength. The K IC values were 
calculated by Eq. (1) using nominal stress. 

The mechanical properties in a single annual ring for 
FEM analyses of sugi, spruce, and akamatsu are shown in 
Tables 2, 3, and 4, respectively. 


Results and discussion 

Figure 9 shows the relations between the location of the 
crack tip in an annual ring and the fracture toughness (A/ IC ) 
values for sugi, spruce, and akamatsu. Experimental values 
of K 1C in Fig. 9 are the same as those in previous reports for 
sugi 1 and spruce. 2 The results for akamatsu are newly ob¬ 
tained here. The variability of K lc values in an annual ring 
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Specific gravity Specific gravity 



Specific gravity Specific gravity 
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Fig. 7. Distributions of the predicted tensile strength in a single annual ring. Filled circles, o tT ; open circles, o m 


Fig. 8. Finite element meshes 
for the single-edge notched type 
tensile specimen of TR crack 
propagation system. Units are 
millimeters 



D E 



was reconfirmed. The FEM analysis results, shown by the 
solid line in Fig. 9, were found to be in good agreement with 
the experimental findings. 

For understanding these K lc variations in an annual ring 
according to the location of the crack tip and the direction 
of crack propagation, we examined stress concentration 
around the crack tip. The stress at the tip is expressed by the 
following equation: 

% = acr (5) 


where the a tjp is a stress at the crack tip in a tangential 
direction obtained by FEM analysis, cr is the nominal stress, 
and a is the stress concentration factor. The results, shown 
in Fig. 10, indicate that a differs depending on the propaga¬ 
tion direction of the crack (pith-side notched, bark-side 
notched), even at the same location in an annual ring. The 
value of a increases with an increase in the relative position 
in an annual ring. The variability of K 1C values derives from 
this tendency of a and the variation of tensile strength 
shown in Fig. 7. 
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Fig. 9. Relations between the K 1C 
value and the location of a crack 
tip in a single annual ring. Left 
Pith-side notched. Right Bark- 
side notched. Circles, K lc by ex¬ 
periment; solid lines, K lc by the 
finite element method 




Relative position of crack tip ( % ) 



Relative position of crack tip ( % ) 



Relative position of crack tip ( % ) 



Relative position of crack tip ( % ) 



Relative position of crack tip ( % ) 


The results of measurements of acoustic emission 
(AE) and fractography of spruce in our previous report 2 
were compared with the results from the present work and 
can be summarized as follows. When the degree of stress 
concentration was small, an AE of a small amplitude was 
gradually generated, and intrawall failure was clearly ob¬ 
served at the fracture surface that was considered to be 
formed before the crack initiation. On the other hand, when 


the degree of stress concentration was large, an AE of large 
amplitude was rapidly generated, and transwall failure was 
clearly observed. Considering the behavior before crack 
initiation, we believe that intrawall failure generated by a 
low level of concentration of stress relaxes stress concentra¬ 
tion at the crack tip and contributes to the increased K lc 
value. 
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Relative position of crack tip ( % ) Relative position of crack tip ( % ) Relative position of crack tip ( % ) 

Fig. 10. Relations between the stress concentration factor (a) and the location of a crack tip in a single annual ring. Filled circles , pith-side 
notched; open circles, bark-side notched 


Conclusions 

The effects of the location of the crack tip in an annual ring 
and the direction of crack propagation on the fracture 
toughness of the TR crack propagation system of coniferous 
wood were analyzed by the finite element method in regard 
to changes in elastic modulus and strength within an annual 
ring. In terms of the variability of K lc values, the experimen¬ 
tal and calculated results had the same tendency. Therefore, 
it was concluded that the variability of K IC is caused by (1) 
the heterogeneity of the elastic modulus and strength within 
an annual ring and (2) the changes in the degree of stress 
concentration at the crack tips, depending on the direction 
of crack propagation. 
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